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SCOPE 


This document records the software which has been developed in response 
to AD004-72-455-02, "Pasture/Range Condition Assessment Project." The 
total software system referred to as the Pasture/Range Condition Assess- 
ment is described in detail. 


SYSTEM OVERVIE W 

The Range processor allows a user to analyze Landsat segment containing 
pasture crop cover. 

The Range processor allows four options; they are: mapping, generating 

statistics on range mask, calculating VINS, and plotting VINS. 

The Range processor can be read on the SPU or the CCA. Plotting is 
allowed on the SPU only. 


V 


1 . RANGE 


Purpose 

This processor is the driver for the Range software. It will allow the user to 
analyze Landsat data containing pasture and rangeland. Analysis includes map- 
ping, generating statistics, calculating vegetative indices, and plotting vege- 
tative indices (VINS). 

Calling Sequence 

Program Range (1.) 

Calling Arguments 

Not applicable 

Files 

Not applicable 
Common Blocks 

The following abbreviations are used in tables throughout this document: 

A = alpha numeric 
I = integer 

I/O = input/output 
L = logical 

R = real 

Common blocks for Range are /RADVAL/, /FILES/, STATS/, and /VINS/. 
/RADVAL/COMMON BLOCK 

Relative 

Name Word Number Usage Type Definition 

IMDATA 1-450 Array I Array for IMDACS image data 


1 


/VI NS/ COMMON BLOCK 


Relative 


Name 

Word Number 

Usage 

Type 

Definition 

AVI 

1 

Process 

I 

Ashburn vegetative index. 

GRN 

2 

Process 

I 

Greeness from Kauth algo- 
rithm. 

ND3 

3 

Process 

I 

ND3 = ^^^^■^*^^V(Ch3+Ch2) 

ND4 

4 

Process 

I 

ND4 = ^^‘^^”^^^V(Ch4+Ch2) 

RAT 

5 

Process 

I 

rat = ^*^*^^/Ch2 

SAVI 

6 

Process 

I 

AVI, corrected for sun angle 

SGRN 

7 

Process 

I 

GRN, corrected for sun angle 

SND3 

8 

Process 

I 

ND3, corrected for sun angle 

SND4 

9 

Process 

I 

ND4, corrected for sun angle 

SRAT 

10 

Process 

I 

RAT, corrected for sun angle 

LAV I 

11 

Process 

I 

AVI, corrected for sun angle 
and satellite 

LGRN 

12 

Process 

I 

GRN, corrected for sun angle 
and satellite 

LND3 

13 

Process 

I 

ND3, corrected for sun angle 
and satellite 

LND4 

14 

Process 

I 

ND4, corrected for sun angle 
and satellite 

LRAT 

15 

Process 

I 

RAT, corrected for sun angle 
and satellite 
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/FILES/COMMON BLOCK 


Page 1 of 2 


Name 

Relative 
Word Number 

Usage 

Jm. 

Definition 

FNAME 

1-30 

Array 

L 

Name of fields file 

GTNAME 

31-60 

Array 

L 

Name of ground truth file 

MNAME 

61-90 

Arra,'' 

L 

Name of mask file 

IMGNAM 

91-120 

Array 

L 

Name of image data file 

/STATS/COMMON BLOCK 




Name 

Relative 
Word Number 

Usage 

lyee 

Definition 

NPRR 

1 

Process 

I 

# of pixels labeled range 
in gt + mask 

NPRN 

2 

Process 

I 

# of pixels labeled range 
- gt, nonrange-mask 

NPNR 

3 

Process 

I 

# of pixels labeled nonrange 
- gt, range-mask. 

NPMN 

4 

Process 

I 

# of pixels labeled nonrange 
ingt + mask 

NPR 

5 

Process 

I 

# of pixels labeled range- 
mask 

NPN 

6 

Process 

I 

# of pixels labeled nonrange- 
mask 

NPMIX 

7-16 

Array 

I 

Array dimensioned (2, 5), used 
in calculating overview statis- 
tics report 

TOTAL 

17 

Process 

I 

Number of pixels in report 

NPIXA 

18 

Process 

I 

Number of pixels that agree 

NPIXD 

19 

Process 

I 

Number of pixels that disagree 

NPIXPA 

20 

Process 

I 

Number of pixels that partially 
agree 
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/STATS/COMMON BLOCK 


Page 2 of 2 


Relative 


Name 

Word Number 

Osage 

IZPe 

Definition 

PROPA 

21 

Process 

R 

Proportion of agreement 

PROPD 

22 

Process 

R 

Proportion of disagreement 

PROPPA 

23 

Process 

R 

Proportion of partial agreement 

PROPC 

24 

Process 

R 

Proportion of pixels committed 

PROPO 

25 

Process 

R 

Proportion of pixels omitted 

PROPNC 

26 

Process 

R 

Proportion of nonrange pixels 
commi tted 

PROPRO 

27 

Process 

R 

Proportion of range omitted 

PROPRC 

28 

Process 

R 

Proportion ‘cf range committed 

PRGTLR 

29 

Process 

R 

Proportion of ground truth 
labeled range 

PRMLR 

30 

Process 

R 

Proportion of mask labeled 
range 

DIFPGM 

31 

Process 

R 

Difference in labeled range 
proportions 
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Subroutines Called 


Name 

Definition 

Reference 

PRSET 

Setup routine, prompt^*- for user input 

1.1 

LBLGT 

Reads digitized ground truth tape and counts 
subpixels, writes map to disk. 

1.2 

GTMAP 

Generates gray scale map file of ground truth 

1.3- 

msmp 

Reads an IMDACS formatted mask file and generates gray 
scale map file. 

1.4 

CMPGTM 

Compares digitized g'^ound truth at subpixel level 
to mask and generates a gray scale map file. 

1.5 

WRTMAP 

Writes a gray scale map printer. 

1.6 

MASKS! 

Compares IMDACS mask file to labeled ground truth and 
computes statistics. 

1.7 

WRSTAT 

Writes the statistics report to the printer 

1 .8 

PRVIN 

Driver for VIN calculations. 

1.9 

PRPLOT 

Driver for VIN plotting. 

1.10 

PLTVIN 

Reads VIN data and plots it on the Cal Comp plotter. 

1.11 


Note that PRPLOT and PLTVIN are used on the SPU. On the CCA, the routines 
are dummy routines. 
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CALLED BY 


RANGE 0.) 

Local Variables 

Name 

Definition 

ACQLST 

Array of acquisition dates. 

AWC 

Not used. 

CARD 

Array for user input. 

CMI 

Crop moisture index. 

CORECT 

Data correction . 

CRYR 

Crop year for ground truth data. 

DATYPE 

Data type. 

DISK 

Data disk. 

GRASS 

Flag for inc] tiding grass crop codes. 

I ERR 

Error flag. 

INBUF 

Input array for data. 

INOPT 

Option. 

LSATS 

Satellite number associated with data on euch acquisition date. 

MOPT 

Option for map type. 

NACQ 

Number cf acquisitions. 

NCOR 

Number of data corrections. 

NDAT 

Number of data types. 

OCTALS 

Octal equivalent of Julian day. 

OTHER 

Flag for CMI data and CMI plots. 

OUTBUF 

Output array for data. 

PIX 

Subpixel counts. 

POPT 

Option for plot type. 

SEGID 

Array used in generating file name. 

SEa^^ 

Segment number 

SLINE 

Line number. 

SOPT 

Option for statistics report. 
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1.1 PRSET 


Purpose . 

This is the setup routine for the Range processor. It queries the user for 
input options. 

Calling Sequence 

Subroutine PRSET (sign, cryr, iopt, thresh, NCC, stype, ccode, NDOT, dotlst, 
segid, mopt, sopt, nacq, acqist, popt, card, grass, disk, ierr). 

Calling Arguments 


Name 

Usage 

Jm. 

Definition » 

Sign 

Input 

I 

Segment number 

Cryr 

Input 

I 

Year of ground truth data 

IOPT 

Input 

I 

Option selector 

THRESH 

Input 

I 

Threshold on subpixel count 


If THRESH = N, then N 

subpixels must be of the crop code 


for the 

pixel to be 1 

abeled as containing this crop. 


IFN = 6, 

the subpixel 

is pure. 

NCC 

I/O 

I 

Number of crops to map.- 

STYPE 

I/O 

I 

Option for 1978 segment, wheat 




or corn. 

CCODE 

Process 

I 

Table of crops to map. 

NDOT 

I/O 

I 

Number of pixels. 

DOTLST 

Array 


List of pixel coordinates. 

SEGID 

Array 

L 

Array used in generating file names. 

MOPT 

Process 

I 

Option for map type. 

SOPT 

Process 

I 

Option for statistics report type. 

NACQ 

Process 

I 

Number of acquisitions. 

ACQLST 

Array 

I 

List of acquisition dates. 

POPT 

Prpcess 

I 

Option for plot type. 

CARD 

Array 

L 

Array for user input. 

GRASS 

I/O 

I 

Indicator to use grass crop codes. 

DISK 

I/O 

L 

Data disk. 

IERR 

I/O 

I 

Error flag. 

Files 




Not applicable 
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Common Blocks 


PRSET uses common block /FILES/, For a description of /FILES/, refer to 
Section 1. 

Functions Called 

Name Definition Reference 

NUMBR Decodes a line of input numbers into an integer array 1.1.1 

DPNUMB . Same as NUMBR, except for INTEGER *4 input values 1.1.2 

NAMEIT Queries for and inputs disk file name 1.1.3 

Called By 

Range (1.) 

Local Variables 

Name Definition 

SEGMT Segment number 

STYPE Indication for 1978 segment for wheat or corn/soybean. 

1. M NUMBR 

Purpose 

This function decodes an array of character numerals into integer 
format. Each number must be separated by a comma. NUMBR return 
the number of integer values in the array. 

Calling sequence 

Function NUMBR (CARD, COL, NUMVEC, VECSIZ) 

Cal 1 ing Arguments 


Name 

Usage 

Type 

Definition 

CARD 

Input 

L 

Array for user input. 

COL 

Input 

I 

Column the first character is in. 

NUMVEC 

Output 

I 

Array of decoded integers. 

VECSIZ 

Output 

I 

The number of integers decoded 
and stored in array NUMVEC. 
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Files 


Not applicable 


Common Blocks 


Not applicable 


Subroutines Called 

Not applicable 


Called By 


PRSET (1.1) 


Local Variables 


Name 

Definition 

CHAR 

Character being decoded 

CRDSIZ 

Set to 80» length of input string 

NUM 

Integer value of CHAR 

SUBTOT 

Integer value of character string, SUBTOT is stored in NUMVEC 

1.1.2 DPNUMB 

. 


This function decodes an array of character numerals into a double preci- 
sion integer array. DPNUMB is a modification of NUMBR. DPNUMB returns 
the number of integer values in the array. 

Calling Sequence 

Function DPNUMB (card, col, numvec, VECSIZ) 

Calling Arguments 


Name 

Usage 

Im. 

Definition 

CARD 

Input 

L 

Array for user input 

COL 

Input 

I 

Column the first character is 

NUMVEC 

Output 


Array of decoded integers 

VECSIZ 

Output 

I 

The number of integers decoded 
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Files 


Not applicable 
Common Blocks 
Not applicable 
Subroutines Called 
Not applicable 
Called By 

PRSET (1.1), PRPLOT (1.10) 

Local Variables 

Name Definition 

CHAR Character being decoded 

CRDSIZ Set to 80 length of input string 

NUM Integer value of CHAR 

SUBTOT Integer value of character stirng, SUBTOT is stored in NUMVEC 

1.1.3 NAMEIT 
Purpose 

This function reads in a disk or tape file name. NAMEIT returns the 
number of characters in the name. 

Calling Sequence 

Function NAMEIT (name, iopt, disk) 

Calling Arguments 


Name 

Usage 

Type 

Definition 

NAME 

Output 

L 

File name 

IOPT 

Input 

I 

Option used to query user 

DISK 

Input 

L 

User UIC 


Files 

Not applicable 
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Common Blocks 


Not applicable 
Subroutines Called 
Not applicable 
Called By 

PRSET (1.1), PRPLOT (1.10) 

Local Variabl es 

Name 

BUFF 
UIC 

1.2 LBLGT 
Purpose 

This routine reads the digitized ground truth file and labels pixels 
according to the map option. For map option 1 or map option 2, subpixel 
crop codes are compared and a count saved for the code that has the great- 
est number. This count is compared to a threshold and a file of counts is 
written. For map option 3, subpixel crop codes are compared to desig- 
nated crop codes and a count is saved. A file of counts is written. 

Calling Sequence 

Subroutine LBLGT (iopt, cryr, gtline, iobuf, thresh, grass, pix) 

Calling Arguments 


Name 

Usage 

Iy-P.e 

Definition 

IOPT 

Input 

I 

Option chosen 

CRYR 

Input 

I 

Crop year 

GTLINE 

Output 

I 

Array for digitized ground truth values 

IOBUF 

Output 

L 

Array of counts 

THRESH 

Input 

I 

Threshold for subpixel count 

GRASS 

Input 

I 

Indicator for including grass crop codes 
with range 

PIX 

ft 

Output 

L 

Array of counts 

n 


Definition 

Buffer for user input 
User UIC, set to [70, 7] 


Files 


File 

Usage 

Record Format 

GTNAME 

Ground truth data 

JSC universal format 

Grayitiap. dat 

Pixel counts 

unformatted records 

SPXCTS. dat 

Pixel counts 

unformatted records 

Common Blocks 




LBLGT uses common block/FILES/. For a description of /FILES/, refer to 
Section 1 . 

Subroutines Called 

Not applicable 

Called By 

Range (l.)» PRVIN (1.9) 
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LQ.c,gJ-.Vflr-lab1 es. 

Name 

CPIX 

FLAG 

INBUF 

INBUFl, INBUF2, INBUF3 

N.CHK 

OFST 

PAS77 

PAS78 

PIXCT 

TOTALS 


Definition 
Subpixel codes 

Flag to indicate if a subpixel has been checked 
Input array for digitized ground truth 
Input arrays 

Number of crop codes to check 
Offset in reading ground truth data 

1977 pasture and grass crop codes 

1978 pasture and grass crop codes 
Counts in output array 
Subpixel counts and codes 


1.3 GTMAP 


Purpose 

This routine creates a gray scale map of ground truth data. 
Calling Sequence 

Subroutine GTMAP (cryr, ccode, iobuf, ncc, stype) 

Calling Arguments 


Name 

Usage 

Jm. 

Definition 


CRYR 

Input 

I 

Crop year 


CCODE 

Input 

I 

Table of crop codes to i 

map 

IOBUF 

Output 

L 

Output array 


NCC 

Input 

I 

Number of crop codes to 

check 

STYPE 

Input 

I 

Crop type for year 1978 

, wheat 




or corn/soy 


Files 





File 

Usage 


Record Format 


MAPOUT.DAT 

Map file 


Unformatted records 


CC0DE77W.DAT 

Crop code symbols 

(13, IX, 13, 2X, Al) 


CC0DE78W.DAT 

Crop code symbols 

(13, IX, 13, 2X, Al) 


CC0DE78C.DAT 

Crop code symbols 

(13, IX, 13, 2X, Al) 


GTNAME 

Digitized ground truth 

JSC universal format 



Common Blocks 
Not applicable 
Subroutines Called 
Not applicable 
Called By 
Range (1.) 
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Local Variables 

Name Definition 

CCTAB Crop codes 

CHECK, CODE Crop code value 

SP Stop pixel number 

ST Start pixel number 

1.4 MASKMP 

Purpose 

This routine generates a gray scale map file of the mask. 
Calling Sequence 


Subroutine MASKMP (outbuf) 
Calling Arguments 


OUTBUF 

Files 


Usage 

Output 


Im. 


Definition 
Output array 


MNAME 


Usage 
Mask file 


MAPOUT.DAT Map file 
Common Blocks 


Record Format 

8 (126A1) header, (58A1) dummy record, 
(126A1), (126A1), (126A1), (32A1) 
data records 

Unformatted records 


MASKMP uses common block /FILES/. For a description of /FILES/, refer to 
Section 1. 

Subroutines Called 

Not applicable 

Called By 

Range 
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Local Variables 


Name Definition 

DUMBl, DUMB2 Arrays used in reading mask file data 
INBUF Array used in reading mask file data 

1.5 CMPGTM 

Purpose 

This routine compares ground truth data to the IMDACS mask file and writes 
a comparison map file. 

Calling Sequence 

Subroutine CMPGTM (lOBUF, GTBUF) 

Calling Arguments 


Name 

Usage 

Tyge 

Definition 

lOBUF 

Output 

L 

Array for output data 

GTBUF 

Input 

L 

Array used in reading counts file 
created by LBLGT 

Files 

File 

Usage 


Record Format 

MNAME 

Mask file 


8 (126A1) header, (58A1) dummy 
record, C126A1), (126A1), (126A1), 
(56A1) data records 

MAPOUT.DAT 

Map file 


Unformatted records 


Common Blocks 

CMPGTM used common block /FILES/. For a description of /FILES/, refer to 
Section 1. 

Subroutines Called 

Not applicable 
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Called By 
Range (1.) 

Local Variables 


NAME 

Definition 



DUMBl . 

Arrays used 

in reading 

input data 

DUMB2 

Arrays used 

in reading 

input data 

INBUF 

Array used 

in reading 

input data 

1.6 WRTMAP 




Purpose 




This routine 

writes a gray 

scale map to the printer. 

Calling Sequence 



Subroutine WRTMAP (Segn, Cryr, lopt, Thresh* lobuf* Cols* Stype* NCG) 

Calling Arguments 



Name 

Usage 

IZfie 

Definition 

SEGN 

Input 

I 

Segment number 

CRYR 

Input 

I 

Crop year 

lOPT 

Input 

I 

Indicator of map type 

THRESH 

Input 

I 

Threshold on subpixel count 

lOBUF 

Output 

L 

Output array for map file. 

COLS 

Output 

I 

Column numbers on map 

STYPE 

. Input 

I 

Indicator for crop type for 19‘ 

NCC 

Not used 



Files 




File 

Usage 


Record Format 

MAPOUT.DAT 

Map file 


Unformatted records 

PRINTER FILE 

Map 


Map is generated on four pages 
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Common Blocks 


WRTMAP used common block /FILES/. For a description of /FILES/, refer to 
Section 1. 

Subroutines Called 

Not applicable 

Called By 

Range (1.) 

Local Variables 


Name 

Definition 

TST 

Row numbers on map 

1.7 MASKS! 



Purpose 

This routine computes mask pixel statistics. 
Calling Sequence 

Subroutine MASKS! (lopt, lobuf, Spixct) 


Calling Arguments 



Name 

Usage 

Type 

Definition 

lOP! 

Input 

I 

Statistics report option 

lOBUF 

Process 

L 

Array used in processing mask file i 

SPIXC! 

Input 

L 

Array used in reading counts file 
created in LBLG! 

Files 

File 

Usage 


Record Format 

MNAME 

Mask file 


8 (126A1) "• header, (56A1) - dummy 
(126A1), (126A1), (126A1), (56A1) - 
data records 

SPIXCTS.DA! 

Counts file 


Unformatted records 
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Common Blocks 


MASKS! uses common blocks /FILES/ and /STATS/. For a description of these 
blocks, refer to Section 1. 


Subroutines Called 


Not applicable 
Called By 
Range (1 . ) 

Local Variables 

Name Definition 


DUMBl, DUMB2 
FLAG 

INBUF 

1.8 WRSTAT 


Arrays used in reading input data 

Indicator that is alternately set to T and F; determines which 
way pixel is thrown if counts value equals 3 
Array used in reading input data. 


Purpose 

This routine writes the statistics report. 


Calling Sequence 

Subroutine WRSTAT (lopt, Segn, Cryr, Grass) 
Calling Arguments 


Name 

Usage 

Jype 

Definition 

lOPT 

Input 

I 

Statistics report option 

SEGN 

Input 

I 

Segment number 

CRYR 

Input 

I 

Crop year 

GRASS 

Input 

I 

Flag to indicate if c'^^ss crop codes are 
included with pasture crop codes 
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Files 



File 

Usage 

Definition 

Printer file 

Statistics report 

Report prints on one page 

Common Blocks 



WRSTAT uses corranon 

blocks /FILES/ and /STATS/. 

For a description of these 


blocks, refer to Section 1. 

Subroutines Called 
Not applicable 
Called By 
Range (1.) 

Local Variables 
Not applicable 
1.9 PRVIN 
Purpose 

This routine is the driver for the option that calculates VINS. 

Calling Sequence 

Subroutine PRVIN(Topt, Segid, nacq, Acqlst, ndot, Dotlst, Inbuf, Outbuf, 
Card, Cryr, Gtline, Thresh, Grass, Octal s, Lsats, Disk, lERR) 

Calling Arguments 


Name 

Usage 

Jm. 

Definition 

lOPT 

Input/Output 

I 

Option for VIN data set 

SEGID 

Input/Output 

I 

Segment identification 

NACQ 

Input/Output 

I 

Number of acquisitions 

ACQLST 

Input 

I 

Acquisition dates 

NDOT 

Input/Output 

I 

Number of pixels 

DOTLST 

Input/Output 

I 

Pixel coordinates list 
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Calling Arguments (Cont^d) 


Name 

Usage 

Jm. 

Definition 

INBUF 

Output 


L 

Array for pixel counts data. 

OUTBUF 

Output 


L 

Array for spectral values. 

CARD 

Process 


L 

Array used in reading character strings. 

CRYR 

Input/Output 

I 

Crop year of ground truth. 

GTLINE 

Output 


L 

Array for digitized ground truth. 

THRESH 

Input 


I 

Threshold for subpixel counts. 

GRASS 

Input 


I 

Indicator for including grass crop codes 
with range. 

OCTALS 

Output 


L 

Octal equivalent of Julian day for each 
acquisition date. 

LSATS 

Output 

• 

L 

Array containing satellite i that acquired 
data for each acquisition date. 

DISK 

Input/Output 

L 

Disk that contains data information. 

lERR 

Output 


I 

Error flag. 

Files 





File 


Usage 


Record Format 

SEGINFO.DAT 


Segment information (80A1) 

FIELDATA.IMD;100 

Spectral 

data 

for Ch 1 (<PIXCT> 13), value of pixel for 

each record is stored in array LL. 

FIELDATA.IMD;200 

Spectral 

data 

for Ch 2 (<PIXCT> 13) 

FIELDATA.IMD;300 

Spectral 

data 

for Ch 3 (<PIXCT> 13) 

FIELDATA.IMD;400 

Spectral 

data 

for Ch 4 (< PIXCT>I3) 


Common Blocks 

PRVIN uses common block /FILES/. For a description of /FILES/, refer to Section 1. 


21 


Subroutines Called 


Name 

Definition 

Reference 

MASKX 

Extracts image data using mask. 

1.9.1 

LBLGT 

Writes a file of subpixel counts 

1.9.2 

GTX 

Extracts image data using ground truth. 

1.9.3 

FIELDX 

Extracts image data using fields file. 

1.9.4 

DOTX 

Extracts image data using dot file. 

1.9.5 

CALVIN 

Calculate VINS. 

1.9.6 

WRTVIN 

Writes VIN file to disk. 

1.9.7 

Called By 



Range (1.) 



Local Variables 



Name 

Definition 


ACQDAT 

Acquisition date. 


ACQUIS 

Array for acquisition dates. 


ANGL 

Sun angle. 


DATFIL 

File name for image data, set in PRVIN. 


lA 

Acquisition counter. 


LINCT 

Number of recoru» in DATFIL. 


LL 

Number of samples in each record in DATFIL. 


LSAT 

Laudsat satellite. 


OPTI 

Option used in LBLGT. 


PIX 

Count data. 


SANGS 

Normalized sun angle array. 


SUNANG 

Sun angle. 


UIC 

User ID set to [70,7] 
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.9.1 MASKX 


Purpose 

This routine extracts pixels from the IMDACS image file as defined by 
an IMDACS mask file. 

Calling Sequence 

Subroutine MASKX (lOBUF, OUTBUF, LINCT, LL, lERR, KCOP) 


Calling Arguments 


Name 

Usage 

Jm. 

Definition 

lOBUF 

Process 

L 

Array for mask data. 

OUTBUF 

Output 

L 

Array for pixel data 
extracted from IMDACS 
image file. 

LINCT 

Output 

I 

Number of records in 
FIELDATA.IMD;N00. 

LL 

Output 

I 

Array containing 
number of samples in 
each record in 
FIELDATA.IMD;N00. 

lERR 

Output 

I 

Error flag. 

KOP 

Input 

I 

Indicator used in 
selecting mask or non 
masked pixels. 

Files 




File 

Usage 

Record Format 

IMGNAM 

IMDACS image file 

See IMDACS documentation for 
format. 

MNAME 

Mask file 

8 (126A1) - header, (56A1) - 
dummy record, 3 (126A1), (56A1) 
data records. 

FIELDATA.IMD-,N00 

Data files 

N = 1 , 2, 3, 4. One file is 
generated for each of four data 


channels. 
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1 ! 


Coirjnon Blocks 

MASKX uses common blocks /RADVAL/ and /FILES/. For a description of these 
blocks, refer to Section 1. 


Subroutines Called 


Name 


Definition 


FILEM, 

Entry AREAD Reads one line in an IMDACS image file 


Called By 
PRVIN (1.9) 
Local Variables 


Name 


Definition 


BBUF 

DUMBl, DUMB2 

I CHAN 

ILINE 

INBUF 

ISTAT 

OFST 

PIXCT 


Array used in reading IMDACS image file. 
Array used in reading mask data file. 
Channel value counter. 

Line counter. 

Array used in reading mask data file. 

Error flag in reading IMDACS image file. 
Offset in reading IMDACS mask file record. 
Counter for samples. 


Reference 


1.9. 1.1 
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1.9. 1.1 FILEM 
Purpose 

FILEM is an assembler language subroutine that reads one line from an 
IMD/ICS formatted image file. FILEM has an ^ntry point named AREAD. 

1.9.2 LBL6T 

Refer to section (1.2) 

1.9.3 GTX 
Purpose 

This routine extracts pixels from an IMDACS image fi’<e as defined by 
the ground truth. 

Calling Sequence 

Subroutine GTX (lOPT, GTBUF, lOBUF, LINCT, LL, lERR). 

Calling Arguments 


Name 

Usage 

Type 

Definition 

lOPT 

Input 

I 

Option used in reading data from 
subpixel counts file. 

6TBUF 

Process 

L 

Array used in reading data from 
subpixel counts file. 

lOBUF 

Output 

L 

Array used in writing FIELDATA. 
IMD;N00 files. 

LINCT 

Output 

I 

Number of records in FIELDATA 
.IMD file. 

LL 

Output 

I 

Array containing number of samples/ 
record. 

lERR 

Output 

I 

Error flag. 
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Files 


File 

Usage 

Record Format 

IMGNAM 

Image data file 

See IMDACS documentation for 
format. 

GTNAME 

Subpixel counts file 

Unformatted records. 

FIELDATA.IMD;N00 

Spectral values 

[ < PIXCT > 13], value of 
PIXCT for each record is 
stored in array LL. 

SPXCTS.DAT 

Not used 


Common Blocks 



GTX uses common blocks /RADVAL/ and /FILES/. 

For a description of these blocks 

refer to Section 1. 



Subroutines Called 



Name 

Definition 

Reference 

FILEM, 
entry AREAD 

Reads one line in an IMDACS image file 1.9. 1.1 

Called By 



PRVIN (1.9) 



Local Variables 



Name 

Definition 


BBUF 

Array used in reading 

IMDACS image file. 

ILINE 

Counter in reading IMDACS image file. 

ISTAT 

Error flag in reading 

IMDACS image file. 

OFST 

Offset in reading IMDACS mask file record. 

OPTION 

Option used in reading subpixel counts file. 

PIXCT 

Counter for samples in 

output data record. 
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1.9. 3.1 FILEM 


Refer to Section (1.9. 1.1). 

1.9.4 FIELDX 
Purpose 

This routine reads IMDACS fields files and extracts field pixels. 
Calling Sequence 

Subroutine FIELDX (lOBUF, LINCT, LL, lERR). 


Calling Arguments 


Name 

Usage 


Definition 

lOBUF 

Output 

L 

Array used in writing 

FIELDATA.IMD;N00 

files. 

LINCT 

Output 

I 

Number of records in 
FIELDATA.IMD file. 

LL 

Output 


Array containing 
number of samples/ 
record 

lERR 

Output 

I 

Error flag. 

Files 

File 

Usage 


Record Format 

IMGNAM 

Image data file 


See IMDACS documenta- 
tion for format. 

FNAME 

Fields data 


Direct access, u'lfor- 
matted records. 

FIELDATA.IMD;N00 

Spectral data 


(< PIXCT > 13), value 
of PIXCT for each reci 
is stored in array LL 
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Common Blocks 


FIELDX uses common blocks /RADVAL/ and /FILES/. For a description of /RADVAL/ 
and /FILES/, refer to Section 1. 

Subroutines Called 


Name 

FLDINT 
FILEM, 
entry AREAD 

Called By 

PRVIN (1.9) 


Definition 

Extracts pixels for field from scanline 
Reads one line in an 
IMDACS image file 


Local Variables 


Name 


Definition 


BBUF 

END 

FLDATA 

FLDHED 

FLDVTX 

FPIX, LPIX 

ICHAN 

ICNT 

IFLD 

ILINE 

INV 

IREC 

ISTAT 

IMAX 

IMIN 

LININT 

MAXLIN 

NINT 

NPIX 

NVER 

OFST 

PIXCT 

VERTXS 


Array used in reading IMDACS image data. 
Highest sample number in scanline. 

Field data array. 

Field data heading. 

Field data vertices. 

Pixel intercepts for given field. 

Counter for data files. 

Vertice counter. 

Field counter. 

Line counter. 

Number of vertices for a given field. 
Associate variable in reading field file. 
Error flag in reading IMDACS image file. 
Maximum line number for a given field. 
Minimum line number for a given field. 

Pixel coordinates for field vertices. 

Number of lines for each field. 

Number of vertices in given field. 

Number of pixels on scanline in given field. 
Number of vertices for each field. 

Offset in storing field vertices. 

Counter in writing scampi es in output array. 
Array containing vertices for one field. 


Reference 

1.9. 4.1 

1.9. 4. 2 
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1.9. 4.1 FLDINT 


Fu rpose 

This routine will return the pixel numbers of those pixels on a given 
line that are contained within the boundaries of a nonreetarigular field. 

Calling Sequence 

Subroutine FLDINT (FIELD, NPTS, FL, YLINE, NSAMP, JJ). 

Calling Arguments 


Name 

Usage 

Jm. 

Definition 

FIELD 

Input 

I 

Nonrectangular field table. 

NPTS 

Input 

L 

Number of vertices. 

FL 

Output 

I 

Array containing pixel intercepts. 

YLINE 

Input 

L 

Scanline number. 

NSAMP 

Output 

1 

Number of pixels in field on scan- 
line. 

JJ 

Output 

L 

Number of pixel pairs in array FL. 


Files 

Not applicable. 
Common Blocks 
Not applicable. 
Subroutines Called 
Not applicable. 
Called By 
FIELDX (1.9.4). 
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Local Variables 


Name 

Definition 


NPTSE 

Number of vertices in field 

minus one. 

NPTSi 

Number of vertices 


RXX 

X intercept. 


RXl, RYl, RX2, RY2 

Field vertices. 


XNMl, YNMl 

Vertice pair before current 

pair. 

XNP2, YNP2 

Vertice pair following currant pair. 

XX 

X integer. 


XI, Yl, X2, Y2 

Field vertices. 


1.9. 4. 2 FILEM 




Refer to Section 1.9. 1.1 
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1.9.5 DOTX 


Purpose 

This routine reads an IMDACS formatted image file and extracts up to ten 
user defined pixels. 


Calling Sequence 


Subroutine DOTX (NDOT, DOTLST, lERR). 


Calling Arguments 

Name 

NDOT 

DOTLST 
I ERR 

Files 

File 

IM6NAM 

FIELDATA.IMD;NjZI0 
Common Blocks 


Usage 


Type Definition 


Input 

Input 

Output 


I Number of pixels to 

extract. 

I Pixel coordinates. 

I Error fag. 


Usage 

IMDACS image file 
Spectral values 


Record Format 

See IMDACS documentation 
for format. 

(< NDOT > 13) 


DOTX uses common blocks /RADVAL/ and /FILES/. For a description of 
these blocks, refer to Section 1. 


Subroutines Called 

Name Definition Reference 

FI LEM, entry AREAD Reads one line in an IMDACS image 1.9. 1.1 

Called By 
PRVIN (1.9) 


Local Variables 

Name Definition 


BBUF 

DOTS 

ILINE 

ISTAT 

PIXEL 

TEMPI, TEMP2 


Image data. 

Ordered pixel pairs. 

Line counts. 

Error flag in reading IMDACS fils. 

Extracted spectral values. 

Used in ordering pixel coordinates in DOTLST. 
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1.9. 5.1 FILEM 

Refer to Section 1.9. 1.1 
1.9.6 CALVIN 
Purpose 

This routine calculates VINS for each pixel. Then the average VIN for 
all pixels is calculated. 

Calling Sequence 

Subroutine CALVIN (lACQ, lOPT, NLINE, SUNCOR, LSAT, LL). 

Calling Arguments 


Name 

Usage 

Jm. 

Definition 

lACQ 

Input 

1 

Acquisition. 

lOPT 

Input 

I 

Data set option. 

NCLINE 

Input 

I 

Number of scanlines 
of data. 

SUNCOR 

Input 

I 

Sun angle correction. 

LSAT 

Input 

L 

Lauu:'r.t number. 

LL 

Input 

I 

Array containing 
number of pixels in 
each scanline. 

File 

Fi 1 e 

Usage 

Record Format 

FIELDATA.IMD;N00 

Spectral values 

(< NPIX > 13), where NPIX 
for each record is stored 
in array LL. 


Common Blocks 

CALVIN uses common blocks /VINS/. For a description of /VINS/, refer 
to Section 1. 

Subroutines Called 

Not applicable 

Called B y 

PRVIN (1.9) 
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Page 1 of 2 


Local Variables 


Name 

Definition 


CHAN 

Spectral data. 


CHANl , 

Arrays used in reading spectral data 

values. 

CHAN2, 

Arrays used in reading spectral data 

values. 

CHAN3, 

Arrays used in reading spectral data 

values. 

CHAN4 

Arrays used in reading spectral data 

values. 

LIA 

Landsat lA to Landsat 2 correction. 


LIB 

Landsat IB to Landsat 2 correction. 


L2C 

Landsat 2C tc Landsat 2 correction. 


L3A 

Landsat 3A to Landsat 2 correction. 


lAVI 

AVI, no corrections. 


IGRN 

Kauth greeness, no corrections. 


IND3 

Normalized difference for Ch 2 and Ch 

3. 

1ND4 

Normalized difference for Ch 2 and Ch 

4. 

I RAT 

Ratio of the Ch 4/Ch 2, 


I LAV I 

lAVI corrected for Landsat 2. 


ILGRN 

IGRN corrected for Landsat 2. 


ILND3 

IND3 corrected for Landsat 2. 


ILND4 

IND4 corrected for Landsat 2. 


ILRAT 

IRAT corrected for Landsat 2. 


ISAVI 

lAVI corrected for sun angle. 


ISGRN 

IGRN corrected for sun angle. 
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Page 2 of 2 


■Local 

Variables 



Name 

Definition 


ISND3 

IND3 corrected for sun angle. 


ISND4 

InD4 corrected for sun angle. 


ISRAT 

IRAT corrected for sun angle. 


PTXCT 

Number of pixels used in VIN calculation. 

RPIX 

PIXCT stored as real value. 


SAC 

Sun angle corrected data. 


TAVI 

Total 

AVI. 


TGRN 

Total 

Kauth greeness. 


TND3 

Total 

normalized difference for Ch 2 and Ch 3. 

TND4 

Total 

normalized difference for Ch 2 and Ch 4. 

TRAT 

Total 

ratio for Ch 4/Ch 2. 


TLAVI 

Total 

AVI corrected for sun angle and Landsat. 

TL6RN 

Total 

IGRN corrected for sun angle and 

Landsat. 

TLND3 

Total 

IND3 corrected for sun angle and 

Landsat. 

TLND4 

Total 

IND4 corrected for sun angle and 

Landsat. 

TLRAT 

Total 

IRAT corrected for sun angle and 

Landsat. 

TSAVI 

Total 

ISAVI corrected for sun angle. 


TSGRN 

Total 

ISGRN corrected for sun angle. 


TSND3 

Total 

ISND3 corrected for sun angle. 


TSND4 

Total 

ISND4 corrected for sun angle. 


TSRAT 

Total 

ISRAT corrected for sun angle. 
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1.9.7 WRTVIN 
Purpose 


This routine writes VIN values as a disk file on EWJ931. 

Calling Sequence 

Subroutine WRTVIN (SEGID, CRYR, lOPT, OCTALS, NACQ, DISK, UIC). 
Calling Arguments 


Name 

Usage 

Im. 

Definition 

SEGID 

Input 

L 

Segment identification. 

CRYR 

Input 

I 

Crop year of ground truth data 

lOPT 

Input 

L 

Data type 

OCTALS 

Input 

L 

Array containing octal equivalent of 
Julian day for each acquisition. 


NACQ 

Input 

I 

Number of acquisitions. 

DISK 

Input 

L 

Data disk label. 

UIC 

Input 

L 

User ID, set as [70,7]. 

Files 




File 

Usage 


Record Format 

VI NAME 

VIN data 


(5F8.3) 


Common Blocks 

WRTVIN uses common block /VINS/. For a description of /VINS/, refer 
to Section 1. 

Subroutines Called 

Not applic:ible. 

Called By 

PRVIN (1.9) 
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Local Variables 


Name Definition 

DATYPE Data type, 

EXTN Extension, set as '.VIN;' 

lACQ Counter for acquisitions. 

VINAME Name of VIN disk file. 

YR Crop year. 

1.10 PRPLOT 

Purpose 

This routine queries user for plotting parameters. 

Calling Sequence 

Subroutine PRPLOT (SEGN, SEGID, CRYR, ACQLST, NACQ, OCTALS, LSATS^ CMI, 
AWC, DATAYPE, NDAT, CORECT, NCOR, CARD, OTHER, lERR). 

Calling Arguments 


Name 

Usage 

Ms. 

Definition 

SEGN 

Input 

I 

Segment number. 

SEGID 

Input 

L 

Segment identification. 

CRYR 

Input 

I 

Crop year of ground truth data. 

ACQLST 

Output 

I 

Acquisition list. 

NACQ 

Output 

I 

Number of acquisitions. 

OCTALS 

Output 

L 

Array containing octal equivalent of 
Julian day for each acquisition. 

LSATS 

Output 

L 

Array containing satellite number 
corresponding to each acquisition 
date. 

CMI 

Output 

R 

Crop moisture index dsta. 

AWC 

Not used 



DATYPE 

Output 

L 

Data type. 

NDAT 

Output 

I 

Number of data types. 

CORECT 

Output 

I 

Correction. 

NCOR 

Output 

I 

Number of corrections. 

CARD 

Process 

L 

Array used in reading character data. 

OTHER 

Output 

L 

Flag that is set if CMI data is read. 

lERR 

Output 

I 

Error flag. 
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Files 

File Usage Record Format 

SEGINFO.DAT I 1 Segment information (80A1) 

PALMER.DAT;! Crop moisture index data (5X, 14. 8X, 12, 

F6.2). 


Common Blocks 


Not applicable 


Subroutines Called 

Name 

Definition 

DPNUMB 

Decodes a string of input numbers into a 
double precision integer array. 

NAMEIT 

Queries for and inputs data file names. 

Called B.V 


Range (1.) 


Local Variables 



Name 

ACQDAT 

ACQUIS 

BDUM 

BEGDAY 

CRD 

DIST 

FLAG 

lACQ 

lAWC 

ICMI 

INPUT 

SEGMT 

VINAME 

WEEK 


Acquisition date. 

Array containing acquisition dates. 

Character data. 

Day 1 for CMI data for years 1974 - 1979. 

Crop reporting district for each acquisition. 
Crop reporting district, 

Flag used in checking plot options. 

Counter in acquisition list. 

Not used. 

Crop moisture index data. 

Character data, 

Segment number as character string. 

VIN data file name. 

Week equivalents of Julian dates. 
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, 12, 37X, 


Reference 

1.10,1 

1 , 10.2 



1.10.1 DPNUMB 

Refer to Section 1.1.2. 

1.10.2 NAMEIT 

Refer to Section 1.1.3. 

1.11 PLTVIN 
Purpose 

This routine plots VIN data on the CALCOMP plotter. 

Calling Sequence 

Subroutine PLTVIN (SEGN, SEGID, CRYR, ACQLST, NACQ, OCTALS, ISATS, 
CMI, AWC» DATYPE, NDAT, CORECT, NCOR, DISK, OTHER, lERR). 


Calling Arguments 


Name 

Usage 

Jm. 

SEGN 

Input 

I 

SEGID 

Input 

L 

CRYR 

Input 

I 

ACQLST 

Input 

I 

NACQ 

Input 

I 

OCTALS 

Input 

L 

LSATS 

Input 

L 

CMI 

Input 

R 

AWC 

Not used 


DATYPE 

Input 

L 

NDAT 

Input 

I 

CORECT 

Input 

I 

NCOR 

Input 

I 

DISK 

Input 

L 

OTHER 

Input 

L 

I ERR 

Output 

I 

Files 

File 

Usage 


VI NAME 

VIN data 

file 


Definition 

Segment number. 

Segment identification. 

Crop year of ground truth data. 

Acquisition list. 

Number of acquisitions. 

Array containing octal equivalent 
of Julian day for each acquisition. 

Array containing satellite number 
corresponding to each acquisition 
date. 

Crop moistuv'e index data. 

Data type(s) 

Number of data types. 

Correction(s) . 

Number of corrections. 

Data disk. 

Flag indicating if CMI data will 
be plotted. 

Error flag. 

Record Format 
(5F8.3). 
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Common Blocks 


Not applicable. 

Subroutines Called 

All subroutines called are Cal Comp plotting routines. See Programming Cal Comp 
Electromechanical Plotters for more information on these routines. 


Name 


Definition 


PLOTS 

FACTOR 

NEWPEN 

PLOT. 

LINE 

SYMBOL 

AXIS 


Initialization routine. 

Scaling factor for graphs, axes and headings. 
Routine for changing pen color 
Positioning routine. 

Writes points or draws lines. 

Writes points or writes heading. 

Draws axis and writes axis scale and label 


Called By 


Range (1.). 


Local Variables 


Name 

Definition 

ATYPE 

Index in choosing data type heading. 

BDUM 

Byte variable. 

DATHED 

Heading for AWC (not used) or CMI data axis. 

DDUM 

Date 

DV 

NACQ + 2 

EXTEN 

Extension, set as '.VIN;' 

FV 

NACQ + 1 

HOT 

Parameter used in placement of plot headings. 

IA,IACQ 

Counter for acquisitions. 

ISYM 

Plotting symbol for points on graph. 

ITYPE 

Counter for data type. 

IVIN 

Counter for plots. 

LDATA 

Vin data, sun angle and Landsat corrected. 

LHED 

Landsat and sun angle heading. 

LSYM 

Graph symbol . 

NOTRAW, 

NOTSUN, 

Flags for data types. 
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Local Variables (Cont'd) 


Name 

NOTSAT 

RDATA 

RHED 

SDATA 

SEGMT 

SHED 

THEAD 

TYPHED 

UIC 

VI NAME 
VINHED 
XDATE 

XDV 

XFV 

XHED 

XPOINT 

YDV 

YFV 

YHED 

YPOINT 

YR 

YRHED 

YVIN 

ZDV 

ZFV 

ZHED 

ZPOINT 

ZVALUE 


Definition 


VIN Data, no corrections. 

Raw data heading. 

VIN data, sun angle corrected. 

Segment number as character String. 

Sun angle heading. 

Data type heading. 

Data correction 

User identification, set as '[70,7]' 

VIN data file name. 

Headings for independent axis. 

Julian dates for acquisitions, initial axis value, and axis 
increment. 

X axis increment. 

X axis initial value. 

Contains 'DATE' as character string. 

X axis data point, initial axis value and axis increment. 

Y axis increment for each graph. 

Y axis initial value for each graph. 

Y heading on CMI plot. 

Y axis data point, initial axis value and axis increment. 
Ground truth year. 

Ground truth year heading. 

Independent axis data. 

CMI axis increment. 

CMI axis initial value. 

CMI axis heading. 

CMI data value, initial axis value and axis increment. 

CMI data. 
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